Assessing the pattern of key factors on women’s empowerment in Bangladesh: Evidence from Bangladesh Demographic and Health Survey, 2007 to 2017–18

Background With half a female population, empowering women can be a key factor in our country’s global advancement. Focusing on household decision-making and attitudes toward wife beating, our study addresses the dearth of research exploring how different socio-economic and demographic factors associated with women’s empowerment evolve over the past decade in Bangladesh (from BDHS 2007 to BDHS 2017–18). Methods Data from four waves of Bangladesh Demographic and Health Survey (BDHS, 2007 to BDHS, 2017–18) were used in this study. We put forth two domains—household decision-making and attitudes toward domestic violence—to assess women’s empowerment. Principal component analysis (PCA) was employed to create women’s empowerment index. To assess the unadjusted association between the selected covariates and women’s empowerment, Pearson Chi-square test and ANOVA F test have been used, while adjusted association has been analyzed through proportional odds model (POM). Results In BDHS 2017–18, women from urban areas experienced ‘high’ empowerment than women in rural areas (56.08% vs. 45.69%). A notable change has been observed in the distribution of women’s empowerment index by education over the survey years. Findings also showed that in all the survey years, division, place of residence, education level, number of living children, media exposure, wealth index, working status, and relationship with household head have been found to have significant association with women’s empowerment index. For instance, women who completed secondary education in 2007, 2011, 2014, and 2017–18, respectively have 14.4%, 31.8%, 24.6%, and 39.6% higher odds of having empowerment compared to those who were uneducated. Further, age at first marriage, spousal age gap, NGO membership etc. emerged as a contributing factor in specific survey years. Conclusion Our study affirmed that, over a ten-year period, women were more likely to protest against physical violence and to participate in various decision-making regarding their personal and social life. Empowerment is notably higher among women in urban residents, those with secondary education, 1–2 children, media exposure, and employment. Policy recommendations should emphasize targeted measures to raise awareness and empower uneducated, unemployed, economically disadvantaged, and physically oppressed women.


Introduction
Women's empowerment is a pivotal focus in global development, recognized as essential for sustainable economic growth and poverty reduction [1].Defined as an interactive process facilitating personal and communal change, empowerment enables individuals to impact institutions affecting their lives and communities [2].Specifically, women's empowerment denotes their full participation in economic, political, and social spheres [2,3], critical for comprehensive social and economic development.
In developing countries, gender inequality poses a significant challenge [4], violating human rights and impeding both women's empowerment and national development.The Global Gender Gap Index, incorporating educational attainment, health, political empowerment, and economic participation, highlights Bangladesh's standing as the 65th country out of 156, indicating persistent gender inequality [5].Notably, Bangladesh excels in women's political empowerment, with women well-represented in national legislatures [6].Despite these strides, challenges persist, as evidenced by the prevalence of gender bias and violence against women, particularly in less-developed countries.
Aligned with international commitments such as the Convention on the Elimination of All Forms of Discrimination against Women (CEDAW) and the Sustainable Development Goals (SDGs), Bangladesh has made notable efforts to eradicate gender-based discrimination and violence [7][8][9].The importance of women's empowerment is underscored by the global Sustainable Development Goal 5, aiming to end gender imbalance and violence against females by 2030 [10].Additionally, the third Millennium Development Goal emphasizes progress toward closing the gender gap and enhancing women's political and economic influence [8].
In Bangladesh, despite these commitments, challenges persist.Gender-based violence, especially within households, remains a concern, affecting approximately 54.2% of married women [11].The sheer size of the female population, constituting 50.46% of the total, accentuates the potential for substantial societal progress through increased women's empowerment [12].
Numerous studies have explored potential variables influencing women's empowerment in Bangladesh and globally [13][14][15][16][17].For instance, studies using data from DHS 2010 and BDHS 2014 proposed indicators like women's participation in household decision-making and opinions concerning wife abuse [13,14].Further, a study emphasized primary components of women's empowerment, including employment status, self-worth, self-confidence, decisionmaking ability, and awareness [15].An evidence-based analysis suggested indicators such as personal freedom, household decision-making, domestic financial decisions, and political independence, exploring the impact of social media, education, community cultural values, women's employment, and household participation rate on women's empowerment [16].Another study focused on women's empowerment through access to health information, finding that urban-dwelling, educated, working, middle-aged women had better decision-making abilities [17].
This study aims to identify factors contributing to women's empowerment in Bangladesh, filling a gap in existing literature.Unlike prior studies, our research uniquely utilizes four waves of BDHS data (2007, 2011, 2014, and 2017-18) to assess temporal changes in covariates impacting women's empowerment.Focusing on household decision-making and attitudes toward wife beating, our study addresses the dearth of research exploring how these factors evolve over the past decade in Bangladesh.Through a comprehensive analysis of evolving factors, we strive to contribute valuable insights into the dynamic landscape of women's empowerment in Bangladesh.

Data
Data from the nationally representative Bangladesh Demographic and Health Survey (BDHS) have been extracted for conducting analysis on women's empowerment.To explore the changes in women's empowerment for one decade, we have used data from four waves of BDHS, collected in 2007, 2011, 2014, and 2017-18.These were nationally representative crosssectional surveys based on a two-stage stratified sample of households.The details of the survey design are described in detail elsewhere [18][19][20][21][22].
In general, in the first stage, a certain number of enumeration areas, EAs (EAs varied for different BDHS) were selected with probability proportional to EA size and with independent selection in each sampling stratum.In the second stage of selection, a fixed number-30 households per cluster-were selected with an equal probability of systematic selection from the newly created household listing.In these surveys, only ever-married women of age 15-49 were interviewed, an exception was found for BDHS 2011 in which ever-married women aged 12-49 were interviewed.But making it analogous to other datasets, we have only taken information about the women who belong to reproductive age 15-49.Missing observations have been eliminated, leaving us with a total of 8931, 16274, 16350, and 18723 observations, respectively.

Ethics approval and consent to participate
The Institutional Review Board of ICF International, Rockville, Maryland, USA (Macro International is now known as ICF International) reviewed and gave their approval to The Demographic and Health Surveys (DHS) Program. 2007Program. , 2011Program. , 2014Program. , and 2017-18 -18 Bangladesh DHS were categorized under that approval.Furthermore, the 2011 and 2017-18 Bangladesh DHS also received approval from another ethical committee: The Bangladesh Medical Research Council.These BDHS were implemented under the authority of the National Institute of Population Research and Training (NIPORT) of the Government of the People's Republic of Bangladesh with financial support from USAID/Bangladesh.Prior to asking, informed consent was obtained from each survey participant.Respondents who refused to consent were not included in the survey.

Outcome variable
Construction of women's empowerment indices.In the conceptual framework (Fig 1 ), we considered a total of nine questions which were classified into two broad dimensions; household decision-making and attitudes toward wife beating [19][20][21][22][23][24].BDHS collected information on household decision-making: (1) their own health care, (2) major household purchases, (3) their child's health care (this information is missing in BDHS 2017-18 data), and (4) visits to their family or relatives.For the above decisions, we made four binary variables (yes/no) where "yes" was defined for the responses of "respondent alone", "respondent and husband/partner" or "respondent and other person", and "no" for otherwise.This index is positively related to women's empowerment and for analysis purposes, these questions were recoded (i.e., 1 for "yes" and 0 for "no").The second index, containing five indicators, was used to gauge women's attitudes toward wife beating.Women were asked if they thought a husband had the right to hit his wife if she (1) refused to have sex with him, (2) left the house without his consent, (3) neglected the children, (4) argued with him, or (5) burned the food.The response "yes" to the above questions means that she believes wife-beating is justified and "no" means she rejects wife-beating for that particular reason.As this index is negatively related to women's empowerment, for the analysis purposes, these questions were recoded oppositely from the first indicator (i.e., 1 for "no" & 0 for "yes").For creating the women's empowerment index (WEI), principal component analysis (PCA) has been used with all nine indicators (4 decisions and 5 reasons) [14,18].Principal components (PCs), with their unique variance features, are normalized linear combinations of random (original) variables.Variance of a principal component is the eigenvalue associated with it.Ratio of k th eigenvalue to the sum of all the eigenvalues represents the proportion of variation explained by the k th component [25].Here, in this study, after PCA, nine PCs (for BDHS 2007 and 2017-18, eight PCs) have been generated among which the first principal component was regarded as the women's empowerment score (WES) since most of the variations in the data were explained by it.The WES was further broken down into 3 equal portions based on quantiles; labeled low, middle, and high for domains below the first, in between the first and second, and above the second quantile, respectively.Finally, our preferred outcome metric is the score index with three ordered categories (i.e., low, medium, and high), where these categories indicate order-wise how empowered a woman is (i.e., low means women have low empowerment) [14,18].

Covariates
On the basis of some previous works of literature, the covariates included in this study are age at first marriage, spousal age gap, respondent's education level, education gap, respondent's current working status, number of living children, religion, number of household members, division, place of residence, media exposure, NGO membership, wealth index, husband's occupation, relationship with household head, and sex of household head [14,[26][27][28][29][30].All of these variables could not be directly extracted from the survey.We had to construct many of these with the help of directly available variables.The definition and their measurements (including categories) of variables are given in Table 1.

Statistical analyses
Simple descriptive, bivariate, and multivariate statistical analyses were performed to achieve different objectives in this study.The univariate analysis was performed to investigate the
Age at first marriage Age at first marriage was categorized into two categories as: below 18 years and 18 or 18+ years [14,27].

Spousal age gap
Spousal age gap was measured by taking the difference between spousal ages [14].

Education gap
The difference between the education levels of the spouses was divided into three categories-no gap, wife with lower education, and wife with higher education [14,27].

Difference between sons and daughters
This variable was made by deducting the total number of daughters from the total number of sons.Negative results indicate a gap in which there are more daughters, positive numbers indicate more sons, and zero means simply no difference between daughter and son [14].

Religion
Islam, Hinduism, Buddhism, and Christianity comprised the four categories of religion in the original dataset.However, for the sake of the study, the religion has been reclassified into two groups: Muslim (Islam) and non-Muslim (Hinduism, Buddhism, and Christianity) [18,19].
Media exposure A woman is exposed to the media if she has access to either of the media sources (Categories: Yes/No) based on the variables of watching TV, listening to the radio, and reading newspapers and magazines [14,26,[28][29][30].
NGO membership Women's NGO membership was divided into two groups: "yes" if they are members of any one of the following organizations: Grameen Bank, BRAC, BRDB, ASHA, Proshika, Mother's Club, or other NGOs; and "no" if they are not [14,28,29].However, because NGO activities are considerably declining in Bangladesh, data on these aspects was missing from the most recent BDHS 2017-18 data.The fact that the government offers more facilities than NGOs could be one of the causes [30].
Husband's occupation Husband's occupation was categorized into seven categories as farmer, labor, service, large business, small business, unemployed and other [14].
Sex of household head Sex of household head was divided into two groups: male and female [19][20][21][22]. https://doi.org/10.1371/journal.pone.0301501.t001 individual frequency percentage of the selected two indicators comprising women's empowerment in the four consecutive BDHS surveys.In the bivariate analysis of women's empowerment with different selected variables, the Pearson Chi-square test and ANOVA F-test have been applied as an attempt to find out the unadjusted association between women's empowerment and selected covariates.The covariates with p-value <0.05 has been considered significantly associated with the outcome variable and obtained significant covariates in at least any one study year have been included in the multivariable regression model.Since we have an ordinal response variable, the ordinal logistic regression model was used to find out the adjusted effects of covariates on women's empowerment.To assess the goodness of fit for the fitted model, the likelihood ratio test (LRT) has also been used [26].All the analyses were performed with the help of two statistical software packages: SPSS version 20, for data sorting, cleaning, etc. and STATA version 14.0, for the analysis purpose.
Ordinal logistic regression model.With qualitative response variables, binary logistic, ordinal logistic, and multinomial logistic regression models can be used for analyzing data.But when analyzing polychotomous data, the ordinal logistic regression model gives a more accurate and efficient estimate of regression coefficients where the response variable acts in an ordinal way with each predictor.In our study, we have response variable (women's empowerment) with three ordered categories (low, medium, and high).Hence, proportional odds model (POM) has been employed in this study.
Suppose a response variable Y with categorical outcomes, denoted by 0, 1, 2,. ..,k, and let x denote a p-dimensional vector of covariates.
If Pr½Y � jjx� ¼ p j is the cumulative logistic distribution function, Then, And the dependence of Y on x for the proportional odds model has the following representation: where, α j is the j th intercept, x ¼ ðx 1 ; x 2 ; . . .; x p Þ0 and b ¼ ðb 1 ; b 2 ; . . .; b p Þ0 is the p×1 vector of unknown regression coefficients corresponding to x [31][32][33][34].

Results
To see how the situation involving women's empowerment improved over the study years, the percentage distribution of variables regarding decision-making and attitudes toward violence has been utilized (Fig 2).In BDHS 2007, 86.8% of women reported having freedom of decision-making, which decreased in 2011 to 82.2%.But in 2014, it again increased about 0.9%, and hit a record high of 88.1% in 2017-18.This implies that women in Bangladesh have acquired greater autonomy in household decision-making.Conversely, the percentage of women who were against violence towards them increased about 28.9% from 2007 to an alltime high percentage of 97.7% in 2011.In the next survey year, it dropped to 71.1%.However, in 2017-18, it increased to 80.4%.While there was a spike in 2011, women's perceptions of the justification of violence have improved across the survey years.Table 2 outlines the percentage distribution of the women's empowerment index (WEI) by selected covariates over the survey years.It reveals that women living in urban areas experienced 'high' empowerment than their counterparts in rural areas in four waves of BDHS.For instance, urban (vs.rural) women who reported a high level of empowerment were 39.56% (vs.29.01%), 36.90%(vs.26.89%), 36.76%(vs.29.12%), and 56.08% (vs.45.69%), respectively, from 2007 survey to 2017-18 survey.Over the survey years, a notable change has been observed in the distribution of WEI by education.To be specific, the highest percentage of women with secondary education getting 'high' empowered was recorded in 2017-18 BDHS (48.05%), followed by 2007 BDHS (33.32%), 2014 BDHS (31.53%), and 2011 BDHS (30.29%).In addition, higher-educated women reported having 'high' empowerment were 3.5 times as having 'low' empowerment (15.81% against 56.05%) in BDHS 2017-18.Moreover, the level of empowerment among women with 1-2 children is higher than that of women with no children.In 2017-18 BDHS, 50.76% of women having 1-2 children experienced high empowerment, which was 34.07%according to BDHS 2007.Similarly, compared to other study years, a greater percentage of women who are exposed to mass media (51.26%) reported having high empowerment in 2017-18 BDHS.However, from the Pearson Chi-square and ANOVA F-test p-values, the study finds significant bivariate association for all the selected covariates with women's empowerment.
A summary of the multivariable analysis results that were obtained by fitting the proportional odds model is compiled in Table 3.The overall goodness of fit of the model for each study survey has been checked through the likelihood ratio test (LRT) based on chi-square distribution.As the p-value for the corresponding LRT was less than 0.001 for each survey model, the LRT ensured that model was well fitted.Multivariable analysis results demonstrate the adjusted odds ratios (OR) for women's empowerment from BDHS 2007 to 2017-18.All the covariates except the sex of the household head were found to have a significant effect on women's empowerment in at least any one study year.The odds of experiencing empowerment differed significantly in all the survey years by division, place of residence, education level, number of living children, media exposure, wealth index, working status, and relationship   with household head.In brief, women living in Sylhet division had lower odds of gaining empowerment as compared to Dhaka dwellers women (OR 2007 : 0.641, OR 2011 : 0.670, OR 2011 : 0.643, OR 2017-18 : 0.765).Women from rural areas had significantly lesser influence (12.6%, 16.2%, 7.9%, and 24.2% lower odds, respectively) in receiving empowerment than those who were from urban areas.There were significant differences in the level of women's empowerment among secondary and higher educated women than uneducated.Women who completed secondary education in 2007, 2011, 2014, and 2017-18, respectively have 14.4%, 31.8%,24.6%, and 39.6% higher odds of having empowerment compared to those who were uneducated.Further, higher educated women in the four BDHS waves, respectively were 2.083, 2.201, 1.794 and 1.935 times as likely to have empowerment compared to those who had no education.Number of living children has been seen as one of the most important covariates in this study with p-value <0.001.The likelihood of having 3-4 children than those couples who had no children from 2007 to 2017-18 BDHS, respectively was 1.978, 3.191, 3.119, and 2.095 times as likely to have empowerment.Similarly, couples with 1-2 children and 4+ children, respectively were more likely to have empowerment than couples with no children.In the present study, women who belong to small size families have a better chance to establish empowerment.Similarly, women from richest families, and are exposed to media have higher odds of being empowered.As expected, the employment status of women in Bangladesh was strongly associated with their empowerment.More precisely, the likelihood of being empowered was 1.271, 1.529, 1.161, and 1.283 times as likely for employed women for the study years, respectively, compared to homemakers.Comparing women who themselves are head of the household with those who have other relationship (i.e., wife, daughter, daughter-in-law, sister etc.) with household head were less likely to have empowerment in household decision-making as well as attitudes toward wife beating.

Discussion
Women's empowerment is a pivotal determinant of sustainable development, attracting global attention due to its impact on societal power dynamics [28,35,36].Numerous studies have explored empowerment through decision-making power or attitudes toward wife beating, and our study contributes by focusing on influential determinants in Bangladesh over a decade [17,18,23,24,[26][27][28][29][30][36][37][38][39][40][41][42][43][44].Our analysis, rooted in previous research, identified key factors associated with women's empowerment.Bivariate analysis highlighted strong associations between several covariates and empowerment, emphasizing the role of marriage age, education, media exposure, employment, and socioeconomic status [14,27,28].Adjusted odds ratios further elucidated these associations, emphasizing their significance in shaping women's empowerment.Urban-rural disparities emerged as a crucial factor, with urban women experiencing greater empowerment, particularly in Dhaka.Education played a pivotal role, corroborating findings from prior studies, showcasing the link between education, decision-making ability, and resistance against intimate partner violence [28,36,[45][46][47].However, the unexpected negative association observed for highly educated wives warrants further investigation.Family dynamics, including the number of children and their gender distribution, revealed nuanced insights.Larger families were associated with reduced empowerment, while the impact of sons versus daughters on empowerment exhibited variations across study years.These findings challenge conventional expectations, suggesting the need for a nuanced understanding of family dynamics in the context of women's empowerment.
Wealth quintiles highlighted socioeconomic disparities, with richer women consistently exhibiting higher odds of empowerment.Religious affiliation, specifically Muslim identity, emerged as a significant factor impacting women's empowerment, echoing the influence of religious beliefs on societal norms and behaviors [38,48,49].Media exposure and employment status corroborated past research, showing positive associations with empowerment [46], while NGO membership emerged as a contributing factor in specific survey years.Spousal occupation also played a role, emphasizing the financial stability's impact on women's empowerment [50].The role of women as household heads showcased a positive association with empowerment, but the low prevalence of female household heads underscored the patriarchal nature of Bangladeshi society.Additionally, regional disparities and variations over the study period highlight the need for targeted interventions based on geographical contexts.

Conclusion
The findings of the study demonstrate that the situation of women's empowerment in Bangladesh has improved over the survey years, as evidenced by the indicators: household decisionmaking and attitudes toward violence.Our study confirmed that women's empowerment is significantly associated with factors such as division, place of residence, education level, number of living children, wealth index, media exposure, working status, and relationship with household head.Additionally, this study also provides evidence that urban dwellers women, secondary educated women, women with 1-2 children, exposed to mass media, employed women reported having more empowerment than the others over the ten-year span of time.Ensuring further steps to open door to women's education, raising awareness among women, and providing more facilities can accelerate women's empowerment rapidly and in a better way.

Limitations
Our study is not beyond limitations.Firstly, this study used data extracted from BDHS having in responses only from women which introduced biases to some extent.Secondly, women's empowerment was defined in our analysis using only two indicators (household decisionmaking and attitudes toward wife beating) but this could be more strongly defined using more indicators like economic decision-making, access to healthcare, physical mobility, the decision on family planning, assets ownership, etc. Due to lack of data, some important covariates were not considered in this study.